Magnetic breakdown in a normal-metal-superconductor 
proximity sandwich 
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We study the magnetic response of a clean normal-metal slab of finite thickness in proximity 
with a bulk superconductor. We determine its free energy and identify two (meta-)stable states, 
a diamagnetic one where the applied field is effectively screened, and a second state, where the 
field penetrates the normal-metal layer. We present a complete characterization of the first order 
transition between the two states which occurs at the breakdown field Hb (T), including its spinodals, 
the jump in the magnetization, and the latent heat. The bistable regime terminates at a critical 
temperature T cr u above which the sharp transition is replaced by a continuous cross-over. We 
compare the theory with experiments on normal-superconducting cylinders. 

PACS 74.50+r, 75.30.Kz 



The superconducting proximity effect in a normal 
metal adjacent to a superconductor has received a revived 
interest in the past decade B. Among the fundamental 
equilibrium problems, the magnetic response of normal- 
metal-superconductor (NS) structures deserves particu- 
lar interest. Experiments have demonstrated the non- 
trivial screening properties of these hybrid structures, ex- 
hibiting a hysteretic magnetic breakdown at finite fields 
as well as a presently unexplained re-entrance in 
the magnetic susceptibility at low temperatures || . The 
investigated samples have typical dimensions comparable 
to the coherence length £jv of the normal metal, attribut- 
ing a key role to the quantum coherence of the electrons 
coupled to the macroscopic phase of the superconductor. 

The self-consistent study of the screening currents in 
a NS sandwich within the framework of the Ginsburg- 
Landau (GL) equation was carried out a long time ago by 
the Orsay group p] . Their work has provided the first un- 
derstanding of the non-linear field phenomena such as the 
magnetic breakdown. However, in the proximity effect, 
the GL equations are at their limit of validity and their 
use is restricted to the dirty limit. The quasi-classical 
Green's function technique allows to describe the clean 
limit using the Eilenberger equations pi and to general- 
ize the dirty limit results using the Usadel equations || . 
Zaikin was the first to derive the magnetic response of 
a normal-metal slab of finite thickness connected to a 
bulk superconductor along these lines jHJ. Most no- 
tably, he found a non-local screening behavior in the 
clean limit linear response which has an appealing simi- 
larity to the one found in superconductors of the Pippard 
type. The applied (static) magnetic field was found to be 
over-screened, the magnetic induction changing sign in- 
side the normal layer. Using numerical methods, Belzig 
et al. 11 1 have investigated the non- linear field regime of 
these equations and found two (meta-)stable mean-field 
solutions in both the clean and the dirty limit. In this 
work, we determine the H-T 'phase' diagram shown in 



Fig. |l| of the normal metal layer in the clean limit, where 
the bistable regime is particularly extended. In thermo- 
dynamic equilibrium, we find a magnetic breakdown at 
Hi, (T) , which is a first order transition separating the 
phase of diamagnetic screening from the phase of mag- 
netic field penetration. 
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FIG. 1. H — T phase diagram of the normal metal slab 
of thickness d. The breakdown field Hb (T) ~ $o/Ajv (T) d 
marks the first order transition between the diamagnetic and 
the field penetration phase (Ajv (T) denotes the penetration 
depth, 3>o the superconducting flux unit). The critical point 
at the intersection of the spinodals H sc (T) ~ $o/d 2 and 
H s h (T) ~ &o/\% (T) separates the first order transition for 
Xn {T) <C d from the continuous cross-over at large tempera- 
ture (Ajv (T) > d). 



Mota et al. [ jl2[|13| ] have recently investigated the mag- 
netic response of metallic cylinders with a superconduct- 
ing core. Their data, which is not described by the results 
of the GL equations, was claimed to be characteristic for 
the ballistic limit (l3). This has motivated us to derive 
the analytic dependence of the clean limit expression for 



1 



the breakdown field Hi, on temperature T and thickness 
d of the normal layer and to compare it to the experi- 
ment. From the free energy of the normal layer, which 
allows us to identify the two (meta-)stable states, we 
determine the spinodals, the thermodynamic breakdown 
field Hi, (T) and find the jumps in magnetization and en- 
tropy at the transition. Furthermore, we obtain a critical 
temperature which marks the upper limit of the bistable 
regime (see Fig. |l|). Finally, we work out the signatures 
of the non-locality in the ballistic regime as they show 
up in the magnetic susceptibility x an d compare them 
with the experimental data. The following discussion is 
divided into four sections, the analysis of the constitutive 
relations (Sec. I), the solution of the magnetostatic prob- 
lem (Sec. II), the determination of the breakdown field 
from the free energy (Sec. Ill), and finally the comparison 
with experiment (Sec. IV). 



I. CONSTITUTIVE RELATIONS 

The quasi-classical Green's function technique provides 
an appropriate description of a metal with nearly spher- 
ical Fermi surface. In a finite magnetic field, the vec- 
tor potential A (x) can be included as a phase factor 
along unperturbed trajectories, provided the dimensions 
of the normal metal are smaller than the Larmor radius 
(tl = PFc/eH, cyclotron radius of an electron travel- 
ing at Fermi velocity). In the ballistic limit, the quasi- 
classical 2x2 matrix Green's function g Uri (x, vp) satisfies 
the Eilenberger equation ||] (e = \e\, % = ks = c = 1), 

- (v f • V)g Wn (x,v_p) = 

[(u>„ + iew F ■ A(x))f 3 + A (x) h,g Un (x, v F )] , (1) 

where the mean-field order parameter A provides an off- 
diagonal potential (fi denote the Pauli matrices, u> n — 
(2n + l)irT are the Matsubara frequencies, [•,•] is the 
commutator). We have excluded elastic scattering pro- 
cesses by assuming T ^> l/r e ; . 

We consider a normal metal slab of thickness d on top 
of a bulk superconductor as shown in the inset of Fig. [y. 
The vector potential A = (0, A (x) , 0) describes a mag- 
netic field B = (0, 0, B (x)) applied parallel to the sur- 
face, which induces screening currents j = (0, j (x) ,0). 
We make the following idealizations in the description of 
the NS sandwich: The superconducting order parameter 
follows a step function A (x) — A8 (—as) (A real), no at- 
tractive interactions being present in the normal layer. 
We assume a perfect NS interface as well as specular re- 
flection at the normal-vacuum boundary. 

In the subsequent analysis we restrict our attention 
to the magnetic response of the normal layer. In the 
proximity effect, the macroscopic coherence of the super- 
conducting condensate induces correlated electron-hole 



pairs in the normal layer through the process of An- 
dreev reflection. The basic process consists of an electron 
traveling forward and a hole traveling backward along a 
quasi-classical trajectory as shown in Fig. || (at discrete 
energies, bound Andreev states are found along these tra- 
jectories). In the presence of a magnetic field, the area 
enclosed by the trajectory (see Fig. ||) is threaded by the 
flux 

$ (a, i?, ip) — <j> A (x) • dx = 2 tani9cos</? J A(x)dx, 

(2) 

which can be expressed through the integral a = 
J Q d A (x) dx times a geometric factor due to the inclina- 
tion of the trajectory (the spherical angles d and <p pa- 
rameterize the direction of the trajectory). The current 
carried along a trajectory depends on the phase factor 
2$ (a, i9, ip) /$ acquired by the propagation of both the 
electron and the hole along the Andreev loop, and we 
arrive at an intrinsic non-local current-field dependence 
j (a). The total current is determined by the sum over 
the currents along the quasi-classical trajectories, 

j (x) = lempF T V (v y Tr [f 3 g Un (x, v F )]) (3) 

7T ' — ' 

UJ n >0 

(the brackets (...) denote the average over the angles 
ip), from which we reproduce the expression first derived 
by Zaikin fuj, 



Ha) = 



w/2 



di) 



tt/2 



dfj (a, i?,<p)) ; 



(4) 



where (a n — 2uj n d/vp cos i9) 



j(0,<p,$(a,0,<p)) 



2epF 



T sin2 & cos 9 ( 5 ) 



w„>0 



A 2 sin 27T$/$ 



^tj„cosha„ + \Jlo\ + A 2 sinha„^ + A 2 cos 2 7r$/$ 

Note that j is independent of x. The induced currents 
for each trajectory depend only on the flux $ modulo 
the superconducting flux quantum $ = irhc/e, reflecting 
gauge invariance. At small fields (a/$ < 1), the current 
response is diamagnetic for all trajectories and the prox- 
imity effect produces screening currents in the normal 
metal. As the field increases to a/$ ~ 1, some of the 
more extended trajectories produce paramagnetic cur- 
rents, since the reduced flux $ £ [— $ D /2, $ /2] they en- 
close becomes negative, and the net diamagnetic current 
response is reduced. As we reach large fields (a/$ S> 1), 
the Andreev levels become mutually dephased due to a 
uniform distribution of the reduced flux. The associated 
currents are randomly dia- or paramagnetic and the net 
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current vanishes. Note that the proximity effect, i.e., the 
existence of the Andreev levels is not destroyed in this 
limit, leading to a finite kinetic energy of the currents 
induced by the magnetic field. 

II. MAGNETOSTATICS 

Owing to the independence of j on x, the Maxwell 
equation —d^A (x) = 4tt? and the constitutive equa- 
tion ([|) combined with the boundary conditions H] 
A (x = 0) = and d x A (x = d) = H can be given a 
formal solution. We arrive at a parabolic dependence 
for A (x) = Hx + 4nj (a) x (d — §) parameterized by 

a = Jq A (x) dx, which in turn is determined through 
the self-consistency condition 

a=— + —j(a)d 3 . (6) 



The total magnetization M. (per unit surface) is defined 
by 

r d 

AttM = / dx (d x A (x) -H) = 2ttj (a) d 2 . (7) 
Jo 

Eq. (||) contains the essential physics of the problem: 
For small fields (a — > 0), the current j —ZH/^nd lin- 
early suppresses the magnetic induction on the geomet- 
ric length scale d to the value B (0) — ► —if/2 at the NS 
boundary. The magnetic induction is thus over-screened 
and assumes an opposite sign at the NS interface. A 
closer look shows that the current is given by the linear 
response expression 

i(a/^«l)^- 47rA ^ (T)d a ; (8) 

which depends on penetration depth Ajv (T) -C d to be 
derived below (Eq. (|l6|)). When inserted back into 
([|), the vector potential is found to be strongly sup- 
pressed to a sa 3H\ 2 N (T) /2, and we obtain a consistent 
diamagnetic solution (i.e., a/$ *C 1) for fields up to 
H < <& /\ 2 N . At large fields, the current vanishes (j — ► 0) 
and the magnetic field penetrates the normal layer. From 
Eq. (^) we find a sw Hd 2 /2, consequently this metallic be- 
havior is expected down to magnetic fields H > $ /d 2 , 
as follows from the condition a/$ Q S> 1 for the Andreev 
levels to be dephased. With $ /d 2 < ^o/Afy the dia- 
magnetic and field penetration solution coexist in the 
regime $ /rf 2 < H < Qo/Xfj. These simple estimates 
for the limits of the bistable regime elucidate the numer- 
ical data of Belzig et al. UM . 



III. THERMODYNAMICS 

In the phase diagram of Fig. |l] the upper and lower 
bounds of the bistable regime found from the above 



mean-field analysis are identified with the spinodals of 
the transition, the super-cooled field H sc ~ & /d 2 and 
the super-heated field H s h ~ <& /\ 2 N (T). In the thermo- 
dynamic equilibrium, a magnetic breakdown occurs at an 
intermediate field, connecting the diamagnetic regime to 
the field penetration regime by a first order transition. 
In the following, we determine this breakdown field and 
the associated entropy and magnetization jump from the 
free energy. 

The energy (per unit surface) of the currents j (x) — 
—SF/5A (x) is obtained via an integration over the non- 
linear current expression, 



F(a) 



= j <y ) da ' 

Jo 



= 4^E 

(uj u cosha„ + y / uj 2 + A 2 sinha„^ + A 2 
(cj n cosha„ + yj lj 2 + A 2 sinh a„J + A 2 cos 2 7r$/$ 



/2 j.tt/2 

d$ / dip sin cost? (9) 
Jo 



log 



F (a) describes the difference in free energy between the 
metal layer under proximity and in the normal state. 
F (a) is a monotonous and strictly positive function, re- 
flecting the absence of condensation energy in the normal 
layer, and expresses the cost of the induced proximity ef- 
fect lying in the kinetic energy of the currents induced 
by the vector potential. The free energy !F(T,J\4) is 
constructed by adding the electro-magnetic field energy 
and subtracting the vacuum field contribution, 



(10) 



We do not include the condensation energy and the ki- 
netic energy of the screening currents in the supercon- 
ductor. The field dependent term of the condensation 
energy might in fact be of the order of the free energy 
in the normal layer and would be expected to produce 
numerical corrections in the results, which are not ac- 
counted for by our idealized choice of the order parame- 
ter A (x) — A9 (—x). The kinetic energy of the screening 
currents ~ H 2 X may be neglected. 

After a Legendre transformation, we obtain the Gibb's 
free energy 



G{T,H) = F (T, M) — MH 



8tt 



(11) 



The field term in Eq. (11) describes the work necessary 
to expel the magnetic field. The extrema of the free en- 
ergy Q with respect to a reproduce the equation of state 
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(0). Fig. || shows the free energy Q (H) as obtained from 
the parameterization of Q and H through a p5| . The 
breakdown field i^j, (T) is determined by the intersection 
of the free energies Q of the two (meta-)stable solutions. 
We note that this procedure is equivalent to the Maxwell 
construction in the magnetization curve M. = —dQ/dH 
of Fig. |. 

In the following, we consider the free energy (^) in the 
two temperature limits T = and Ta <C T < A and 
obtain (Ta — Vf/2ttcI), 



F t =o (a) 



Ft»t a (a) 



Pf 



An 3 dm 



tt/2 



d-d 



tt/2 



dip sin $ cos 2 $ 



{arctan[tan7r$/$ ]} , 
2 -7 2 (T,A)/ dd 



(12) 



±P F T 2 



tt/2 



exp 



2T 



T 4 cos 



sin 2 7r<i>/$o 



dip sin cos d 



(13) 



The finite value of the superconducting gap A is ac- 
counted for by the dimcnsionless parameter 



1 (T,A) = A/[JA 2 + (nT) +* T < L (I 4 ) 



0.04 




Hd/$ 

FIG. 2. Magnetization M (H) and free energy Q (T, H) at 
a temperature T = 1.5 Ta [Ta = VF/2nd). The representa- 
tion is universal in the thickness d. The small and large field 
branches represent (meta-) stable solutions describing the dia- 
magnetic and field penetration phases, which overlap in the 
field interval H S c < H < H s h- The Maxwell construction 
determines the first order transition between the phases at 
the breakdown field Hi, (dashed line). Inset: Cross-section 
of the normal-metal slab in contact with the bulk supercon- 
ductor. The quasi-classical electron-hole trajectory at angles 
■d, ip = encloses a flux that enters as a phase factor in the 
wavefunction. 

The free energies of the two (meta-)stable stable states 
can be approximated by their asymptotic forms in the 



limits a — > and a — > oo. In the diamagnetic regime, the 
expansion in a/$ up to quadratic order of ( |l2| ) and jl^ ) 
provide the result 



F(a) 



8tt\ 2 n (T)d' 



(15) 



Eq. (g_5() is valid both in the low and high temperature 
limits using the penetration depth 1 / A 2 ^ (0) = 1/ A 
(47rne 2 /m) at T = and 



(T) 



X 2 



y 2 (T,A)^i e - 2 ^ 



(16) 



for T ^ Ta- Note that the derivative j = —dF/da ap- 
plied to Eq. ( |l5| ) produces the linear response constitutive 
relation of Eq. (|h The Gibb's free energy follows from 
Eqs. jl0| ) and (15]), using the solution of the Maxwell 
equations, 



£(a< *□) 



32tt 



-H 2 d. 



(17) 



Eq. ( |l7j ) is dominated by the magnetization work neces- 
sary to expel the field, which is parametrically larger (by 
(d/Ajv (T)) 2 ) than the kinetic energy of the currents. 

In the field penetration regime we approximate the free 
energy by its asymptotic value at a — > oo. In this limit 
we replace the strongly oscillating functions of <E> in (|l2] ) 
and (|l^ ) by their average value ((arctantan $) 2 ) = 7r 2 /12 
and (sin 2 3>) = 1/2 and obtain 



Qt=o (a > *□) 



Gt»t a (a > $o 



1 $ 2 
384tt \ 2 N d' 

7 2 (T,A) 



16tt 3 



(T>- 

° f 

A 2 ^ 



-2T/7U 



(18) 



The magnetization energy vanishes in this limit. The 
corrections to the free energy ([l8]) are of relative order 
(*o/a) 2 . 

The magnetic breakdown field Hb (T) is determined by 
the intersection of the two asymptotics of the free energy 
g given by Eqs. @ and @, 



# & (T = 0) 



1 



6 Awii 
V2 

7T 



(19) 



4>„ 



ff b (T » T A ) « —7 (T, A) ^ e - d /««( T ). (20) 



We note three important features of this result: The 
temperature dependence is a simple exponential with the 
exponent d/£ (T) = T/Ta, where £iV (T) = v F /2nT de- 
notes the normal metal coherence length. The amplitude 
of the breakdown field scales inversely proportional to 
the thickness of the normal layer, Hi, ~ 1/d. In the limit 
T — > the magnetic breakdown field saturates to a value 
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which is suppressed by the universal factor « 0.37 

as compared to the extrapolation of the high temperature 
result. 

We arrive at the H ~ T phase diagram shown in Fig. 
[IJ. The first order transition between the diamagnetic 
and the field penetration regime takes place between the 
spinodals H sc ~ $ /d 2 < H b (T) < H sh ~ <I> /\ N (T) 2 
which delimit the (meta-)stable regime. Their intersec- 
tion marks the critical temperature 



Tent ~ T A log (d/X N ) . 



(21) 



where Xn (T cr it) ~ d. Below T cr n the penetration depth 
is small, \n {T cr it) < d, and we observe a first or- 
der transition. Above the critical point T cr it, where 
A at (T cr it) > d, a continuous and reversible cross-over 
between the diamagnetic and field penetration regime is 
expected. We note that this distinction is similar to the 
one between Type I and Type II superconductors with 
respect to the penetration depth A and the superconduct- 
ing coherence length £. 

The latent heat (at T ^> Ta) of the transition follows 
from Eqs. @, @, and M) using S = -dQ/dT, 



TAS : 



3 T 



H%{T)d, 



16ttT a 

and is related to the magnetization jump 



4ttAM w -H b (T) d 



(22) 



(23) 



via the Clausius-Clapeyron equation. 

In the derivation of the breakdown field we have used 
the asymptotic expansions of the free energies in a/$ 
and $ /a, respectively. Their quality at the transition 
point is determined by the range of overlap between the 
diamagnetic and the field penetration regimes in Fig. |^, 
which is governed by the parameter A at (T) jd. In the 
diamagnetic phase, the corrections are of the order of 
(a/$ ) 2 ~ (H b \%(T)/<S> ) 2 ~ (X 2 N (T)/d) 2 , and simi- 
larly in the field penetration regime. The expansion thus 
breaks down at A at (T) ~ d, which is the critical point 
of the transition line. We note that the total magnetiza- 
tion changes from its diamagnetic value M. ~ H b d to the 
strongly suppressed value M ~ H b d(\N (T) / d) 2 at the 
transition, reflecting its strong first order character. 



IV. EXPERIMENT 



Mota et al. |12|,|13j have measured the breakdown field 
in Ag-Nb cylinders. The clean limit theory valid for 
T ^> VF/lei may be used provided that l e i 3> d, which 
is claimed to be fulfilled in the experiment. In our com- 
parison we neglect the influence of diffusive boundary 
scattering or any potential barrier at the NS interface, 



and ignore the difference in geometry, cylindrical for the 
sample and planar in the theoretical model. 
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FIG. 3. Breakdown field Ht (T) and linear susceptibility 
47i~x from theory and experiment (the analysis applies to a 
Ag-Nb sample of thickness d = 5.5^im). Theory: results of 
Eqs. (^) and (^i|) shown as solid lines; Hb (T) is rescaled to 
fit the zero temperature value Hi, (0) (horizontal line) to the 
experiment. Experiment: data shown as solid dots, the dot- 
ted line is a guide to the eye C3. Note that the logarithmic 
slope of the breakdown field is reproduced precisely (dashed 
line), while the one of the susceptibility is much smaller than 
expected. 

In Fig. H we show the two data sets for the breakdown 
field data obtained on heating and cooling a sample of 
thickness d = 5.5/im exhibiting hysteresis (the theoreti- 
cal values of the super-cooled and super- heated fields H sc 
and H s h are not reached in the experiments). The data 
saturates at low temperatures, in qualitative agreement 
with our theoretical analysis. Given the electron density 
in Ag, n = 5.8 • 10 22 (=> X N = 2.2 • 10~ 6 ) and d = 5.5^m, 
the breakdown field is determined by Eq. ( |l9|) and d2C|). 
Due to the idealization of our model, which assumes a 
step function for the order parameter, as well as the dif- 
ference between the planar and cylindrical geometry, we 
expect a numerical factor correcting the amplitude of H b . 
Making use of the scaling factor w 0.56 in Eqs. jl9| ) 
and (p0|), we calibrate the theoretical result to fit the 
zero temperature value H b (0), as shown in Fig. |3| The 
theoretical prediction for the high temperature behavior 
then follows from Eq. ( p0| ) and is shown as a solid line 
in Fig. ||. Most importantly, Eq. ( p0| ) accurately repro- 
duces the logarithmic slope —\/Ta of the experimental 
data, thus correctly tracing the signature of the Andreev 
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levels. The amplitude of H\, (T) deviates from the data 
by the constant ratio w 0.64, which can be attributed to 
the presence of a barrier at the NS interface, see below. 

An important further agreement between theory and 
experiment is found in the scaling of the breakdown field 
with sample thickness d, which was reported to be oc 1 /d, 
in accordance with Eq. (20|) (the experimental study in- 
volved 10 samples |l3|,[l6] with thicknesses ranging from 
d = 2.9/im to d = 28/im). Similarly, the critical tempera- 
ture determined in the experiment exhibits the same 
scaling oc 1/d, in agreement with Eq. ( fjl| ) (Ta oc 1/d). 

The only result in the literature on the breakdown field 
H was derived from the GL equations in the dirty limit 
Li <C £d = \J vpl e i/6TvT < d, with the coherence length 
£d limited by Ad (d) <C £d < d (Ad {%) is a space and 
temperature dependent penetration depth, see Ref. 0). 
The breakdown field 



H D (T) « 1.9 



4>c 



Ad (0)£d 



exp (~d/£ D ) 



(24) 



exhibits a simple exponential dependence ond/^D oc vT, 
the amplitude being temperature independent. Further- 
more, no dependence of the amplitude on the thickness 
is present. Clearly, the experimental data deviates sig- 
nificantly from the predictions made by the GL theory. 

The good agreement between the clean limit theory 
and experiment for the breakdown field does not triv- 
ially generalize to other physical quantities, however. In 
particular, the temperature dependence of the linear sus- 
ceptibility x = ■M./H exhibits distinct features due to 
the non-locality of the constitutive relation j (a) which 
are not observed in the experiment. From Eqs. (]?]) and 
(0) we obtain the susceptibility 



4ttx 



A-kM 
H 



3 

'Z r 



i 



3X% [T)/d?- 



(25) 



which exhibits a temperature dependence much like a 
Fermi-Dirac distribution: Attx decays exponentially oc 
\/X 2 N (T) at large temperatures, twice as fast as the 
breakdown field. 



imal value at 2>X 2 N 



The susceptibility takes half its max- 
Ti/ 2 ) /d 2 ~ 1, which roughly coin- 
cides with the critical temperature T cr u (see Fig. |^). 
Below the critical point, the susceptibility saturates as 
the penetration depth decreases below the sample thick- 
ness (Aat (T) < d). Due to the non-locality, the penetra- 
tion depth drops out of the expression for Attx ~ — 3/4 
and we are in the regime of over-screening. The log- 
arithmic derivative at T — T1/2 is predicted to be 
X' (Ti /2 ) fx (T1/2) = l/T A . In Fig. | we show the lin- 
ear susceptibility according to the clean limit predictions 
( p5| ) (there is no fitting parameter). The experimental 
data fails to show the typical saturation of the suscep- 
tibility expected below the critical temperature. At low 
temperature the experimental value clearly exceeds the 
maximal diamagnetic value —3/4 found in the clean limit 



(note that we do not consider the anomalous re-entrance 
effect of these samples at very low temperatures here). 
Most strikingly, the decay at large temperature is slower 
than the decay of the breakdown field, while Eq. ( p5| ) pre- 
dicts a decay with twice the logarithmic slope, see Fig. 0. 
Thus the magnetic behavior of the quasi-ballistic samples 
deviates from the clean limit theory, indicating that the 
elastic mean free path is of order of the thickness of the 
sample. We find that the susceptibility emerges as a very 
sensitive indicator of the non-locality of the constitutive 
relations. 

Let us address the question whether the consideration 
of a insulating barrier at the NS interface may lift the 
discrepancy. The consequences of a finite reflectivity 
at the NS interface on the linear current response has 
been analyzed by Higashitani and Nagai jL7|. Their re- 
sults allow for the reflection coefficient R to be included 
in the penetration depth Aat (T) by redefining the fac- 
tor 1R {T A « T < A) = (1 - R) j (1 + R), in Eq. ©; 
Aat (0) = Aat remains unchanged jpjfl- Inserting the 
modified penetration depth into Eq. (|25|) we obtain the 
linear susceptibility. The additional factor 7 does not 
change the characteristic shape of the susceptibility (sat- 
uration, logarithmic slope at T\i%, exponential decay), 
but only lowers the position of the half-value of x to 
T 1/2 w log [d (1 -R)/X N (l + R)]. Thus the finite reflec- 
tion does not remedy the qualitative discrepancy between 
theory and experiment, which may have to be attributed 
to disorder or diffusive boundary scattering. 

Considering the structure of the equations we may ex- 
pect the dependence on the reflection R to enter in a sim- 
ilar fashion into the breakdown field H\, (T) , although we 
note that this has not been shown rigorously. Eq. ( |25| ) 
inserted in Eq. ( p0| ) gives the high temperature behavior, 
while the zero temperature result of Eq. (^|) remains 
unchanged. We fit the breakdown field data by using 
first an overall scaling factor needed to adjust Hi, (0) and 
secondly, a finite reflectivity, which only enters at high 
temperatures. The fit of the high temperature behavior 
provides us with an estimate of the reflectivity R « 0.21, 
and is represented by the dashed line in Fig. |3j. 

In conclusion, we have calculated the clean limit ex- 
pression for the breakdown field separating the diamag- 
netic phase and the field penetration phase by a first or- 
der transition. We have determined the spinodals, the 
critical temperature as well as the latent heat of the 
transition. In comparison with the experimental data 
on quasi-ballistic samples, we have found good agree- 
ment with respect to the dependences on temperature 
and thickness of Hf, (T, d) and T cr n (d). The inclusion of 
a finite reflection at the NS interface permits an accurate 
fit of the breakdown field and gives an estimate for the 
quality of the NS interface. However, with regard to the 
linear susceptibility, the experiments disagree with the 
clean limit theory, showing the need to include additional 
scattering processes. The susceptibility thus emerges as 
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a quantity which is very sensitive to the non-locality of 
the constitutive relations. 
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